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Abstract 

In the heavy quark limit, both vector mesons in the charmless B — > VV decay should have zero 
helicity and the corresponding amplitude is proportional to the form factor difference {A\ — A2). 
The first observed charmless B — » VV mode, B — ► <\>K* , indicates that the form factors Ai(q 2 ) 
and A2(q 2 ) cannot be very similar at low q 2 as shown in some form- factor models. The approach 
of QCD-improved factorization implies that the nonfactorizable correction to each partial-wave or 
helicity amplitude is not the same; the effective parameters a% vary for different helicity amplitudes. 
The leading-twist nonfactorizable corrections to the transversely polarized amplitudes vanish in the 
chiral limit and hence it is necessary to take into account twist-3 distribution amplitudes of the 
vector meson in order to have renormalization scale and scheme independent predictions. Branching 
ratios of B — > VV decays are calculated in two different models for form factors, and the predicted 
decay rates are different by a factor of 1.5 ~ 2. Owing to the absence of (S — P)(S + P) penguin 
contributions to the Vk-emission amplitudes, tree-dominated decays tend to have larger branching 
ratios than the penguin-dominated ones. 
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1. It is known that the decay amplitude of a B meson into two vector mesons is governed 
by three unknown form factors Ai(q 2 ), A 2 (q 2 ) and V(q 2 ) in the factorization approach. It 
has been pointed out in ||] that the charmless B — > VV rates are very sensitive to the form- 
factor ratio AijA\. This form-factor ratio is almost equal to unity in the Bauer-Stech-Wirbel 
(BSW) model 0, but it is less than unity in the light-cone sum rule (LCSR) analysis for 
form factors ||. In general, the branching ratios of B — > VV predicted by the LCSR are 
always larger than that by the BSW model by a factor of 1.6 ~ 2 |T|. This is understandable 
because in the heavy quark limit, both vector mesons in the charmless B — > VV decay 
should have zero helicity and the corresponding amplitude is proportional to the form factor 
difference (Ax — A 2 ). These two form factors are identical at q 2 = in the BSW model. We 
shall see that the first observed charmless B — > VV mode, B — > <f)K*, recently measured 
by CLEO §, BELLE [§ and BABAR [§, clearly favors the LCSR over the BSW model for 
B — V transition form factors. 

In the present paper we will embark on a study of B — > VV decays in the approach 
of QCD-improved factorization which enables us to compute nonfactorizable corrections in 
the heavy quark limit. In the so-called generalized factorization approach, it is assumed 
that nonfactorizable effects contribute to all partial-wave or helicity amplitudes in the same 
weight. We shall see that it is not the case in QCD factorization. Moreover, we will show 
that the leading-twist nonfactorizable corrections to the transversely polarized amplitudes 
vanish in the chiral limit. Hence it is necessary to go beyond the leading-twist approximation 
for transversely polarized states. 

2. In general the B — > VV amplitude consists of three independent Lorentz scalars: 

A[B(p) -> V 1 (e 1 ,p 1 )V 2 (e 2 ,P2)] oc e^ef^g^ + bp^p„ + ice^apPiPz), (!) 

where the coefficient c corresponds to the p-wave amplitude, and a, b to the mixture of s- 
and o?-wave amplitudes. Three helicity amplitudes can be constructed as 

H 00 = — — \{m 2 B - m\ - m\) a + 2m 2 B p 2 c b\ , 

2mim 2 L J 

H±± = a =F m B p c c, (2) 

where p c is the cm. momentum of the vector meson in the B rest frame and mi (mj) is the 

mass of the vector meson V\ (Vz). For H to occur quark spin in the emitted vector meson 

V 2 has to be flipped. Therefore, the amplitude H is suppressed by a factor of m 2 /m B 0. 

The H ++ amplitude is subject to a further chirality suppression of order mi/mg. In general, 
it is thus expected that |i?oo| 2 > \H | 2 > \H ++ \ 2 . The total decay rate is given by 

T{B - VM) = ^V(|#oo| 2 + \H ++ \ 2 + |#__| 2 ). (3) 
87im B 

In terms of the decay constant and form factors defined by @: 
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(V(p,e)\V,\0) = f v m v e;, 

(v(p\e)\v,\p( P )) = — \ — e^yyv^ 2 ) 

VTip + my 

(V( P ',e)\A^\P(p))= i {(m P + m v )e;A 1 (q 2 ) 



e ■ p 



mp + my 



{p + p')^A 2 {q 2 



-2m\ 



■ p 



[A 3 (q 2 )-A (q 2 )}}, 



where q = p — p', A 3 (0) = Aq(Q), and 



A 3 (q 2 ) = mp + mV A 1 {g 2 )- mp ~ mV -i <- 



2my ' 2my 

one has the factorizable B — > V\V 2 amplitude: 



■A 2 (g 2 



-ifv 2 m 2 



(el-e* 2 )(m B + mi )A^(m 2 2 ) 



2A^\ 



fci • P fl )(e 2 • P s h 7 + »W £ 2 £ i P Pi 

(mfl + mi) 



2y Wl (ml) 
K + mi) 



(4) 



(5) 



(6) 



Take the decay I? — > as an example. In the naive factorization approach for hadronic 
weak decays, the decay amplitude of B~ — > <pK*~ reads (in units of Gpj\p2) 



a 3 + a 4 + a 5 



--(a 7 + a 9 + a 10 



xf^,« + (fl4 + aio)x r^ 

2(a 6 + a 8 )(A— 0|s(l + 75 )u|0) (0|u(l - ^)b\B~) 



(7) 



where a 2i = c 2i + ^-c 2i _i, a 2i _i = c 2i _i + ^-c 2i . Neglecting the annihilation contributions 
from the last two terms in Eq. (|7|), we obtain 



H 



oo 



~a{<t>K*)U 



fn B + m K - 



A 2 ( m l) 



[m B + m K *)A? K \ml) T J™*^ V BK \m\) 



m B + tuk* 

where a((f)K*) = a 3 + a 4 + a 5 — Ma-j + a 9 + ai ). In the heavy quark limit, it is clear that 



H 



00 



[A* K *(0)-A»«*(0)), 



BK* 



2mjc* 

H ±± = ^K*)Um^m B [Af K \0) T V BK \0)]. 



(9) 
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In the so-called generalized factorization, nonfactorizable effects are parametrized in 
terms of N* , the effective number of colors. This amounts to assuming that nonfactor- 
izable corrections weight in the same way to all partial-wave or helicity amplitudes. For 
example, the coefficient a(cj)K*) appearing in Eq. (Q) is postulated to be the same for S, P 
and D (or ifoo and H±±) amplitudes after including nonfactorizable contributions. Clearly 
there is no any known physical argument for justifying this assumption. 

Fortunately, the QCD-improved factorization approach advocated recently in [[| allows 
us to compute the nonfactorizable corrections in the heavy quark limit since only hard inter- 
actions between the {BV\) system and V2 survive in the — > 00 limit. Naive factorization 
is recovered in the heavy quark limit and to the zeroth order of QCD corrections. In this ap- 
proach, the light-cone distribution amplitudes (LCDAs) play an essential role. The LCDAs 
of the light vector meson of interest are given by |T0|J9|] 



(V(p,e)\q a (yWp(x)\0) 



fvm v 



du e i< ^ u P-y +u P- x ) 



\®Y\{u)+fJl\u) 



+ 



+ 



^ fv 



fv m Ql + 771 



''12 



*W<i a) («) 



a/3 



due i{up - y+Up - x \£*J) af} <&\{u) 



(10) 



where z = y — x with z 2 = 0, eff (ej_) is the polarization vector of a longitudinally (trans- 
versely) polarized vector meson, u is the light-cone momentum fraction of the quark q in 
the vector meson, u — 1 — u, fv and fy are vector and tensor decay constants, respectively, 
but the latter is scale dependent. To a good approximation one has ejj 1 = p v /m v for a light 
vector meson. It follows that p ■ e± = 0. It is easily seen from Eq. (|1|) that $11 contributes 
to S and D amplitudes, while $^ to P and S waves. In Eq. (|10D , $||(u) and $j_(m) are 
twist-2 DAs, while and are twist-3 ones. Note that twist-3 longitudinally polarized 
distribution amplitudes and [W\ are not shown in Eq. (P~D|) . The reason for keeping 
the twist-3 transversely polarized distribution amplitudes g^'^ rather than the longitudinal 

(s t) 

ones h\, ' will become clear shortly. 



3. We will now study charmless B — > VV decays within the framework of QCD-improved 
factorization. The power corrections such as the annihilation diagrams can be neglected in 
the heavy quark limit. It turns out that the twist-2 DA contributions to the vertex 

corrections and hard spectator interactions vanish in the chiral limit.Q Hence, we will work 



*For the color-suppressed mode B — * J/ipK*, the twist-2 transverse polarized distribution am- 
plitude is on the same footing as the longitudinal one since J/ip is not massless in heavy quark 
limit. 
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to the leading-twist approximation for longitudinally polarized states and to the twist-3 
level for the case of transverse polarization. As discussed before, the effective parameters a, 
entering into the helicity amplitudes H 0Q and H±± are not the same; they are given by 

+(<* + <*) £ G*( Sj ) + |(c 8 + do) 2 CiG'ts.) + ^[ejG'ts,) - iCfe)] + c„G; ^ 
«i = * + £ + ££*(-f*-12), (11) 



a 3 



ft. 



a 6 - c 6 + — , 



C5_ 



a 7 



-+t + Sf C8( - Fk - 12) -^ c * 



ft C 7 

a 8 =c 8 + -, 
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ft 



cio a s C F , a h 

iV c 47T iV c 97T 

c 9 a s C F a 



ai °- Cl0 + iv; + 4^iv: C9F "o^' 

where CV = (iV^ — 1)/(2N C ), Sj = mf/m^, A g = V q bV* ql , q' = d,s and the superscript h 
denotes the polarization of the vector mesons: /i = for helicity 00 states, and h = ± for 
helicity ±± states. 

There are QCD penguin-type diagrams induced by the 4-quark operators O, for i = 
1, 3, 4, 6, 8, 9, 10. The corrections are described by the penguin-loop function G h (s) given by 

G°(s) = - - - In — + 4 / dx&fAx) f 1 duu{\ - u)\n\s - xu(l - u)l 

3 3 777,6 Jo JO 

G ± (s) = ^-^ln-^+4 f 1 dxg ( l\x) duu{\ - u) ln[s - xu(l - u)l (12) 

3 3 771;, JO JO 

T l f'dx f 1 duu(l — u)l — 2 In — — h ln[a - - it)] + x < 1 ~ u "> \ 

2 Jo x Jo I m 6 1 v ;j s _ OT (i _ M ) J 

In Eq. (O) we have also included the leading electroweak penguin-type diagrams induced 



by the operators 0\ and 2 [11 



£? = (^ W + ^ W )( (S + ^). (13) 
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The dipole operator O g will give a tree-level contribution proportional to 

G> . 2 r dx m, Gt= . 2 f dx m + h lT hf dx m. (14) 

Jo x a Jo x 2 2 Jo x 

In Eq. (|TTD, the vertex correction in the naive dimensional regularization (NDR) scheme 
for 75 is given by 

ir k = -121n-^--18 + tf + (15) 

where the hard scattering function // arises from vertex corrections and /// from the hard 
spectator interactions with a hard gluon exchange between the emitted vector meson and 
the spectator quark of the B meson. Note that the twist-2 transversely polarized distribu- 
tion amplitude <3>^ does not contribute to F: it does not give rise to the scale and scheme 
dependent terms — 121n(/x/mft) — 18 and the hard scattering kernels ff and ffj are propor- 
tional to the light quark mass and hence can be neglected. Therefore, the parameters af 
at the twist-2 level are not renormalization scale and scheme independent. Consequently, it 
is necessary to take into account the twist-3 effects for transversely polarized vector meson 
states. This is why we keep twist-3 DAs ^j* in Eq. (|T0|). As for the helicity zero case, 
it is dominated by the leading-twist one and hence the twist-3 DAs /i| , which are power 
suppressed by order of my/ms, are not considered there. It should be stressed that ag 8 
do depend on the choice of the renormalization scale and scheme. Their scale and scheme 
dependence is compensated by the corresponding (S — P)(S + P) hadronic matrix elements. 
An explicit calculation for fj yields 

" x . , / 1 - 2x 



/° = J dx®\{x) (z- ^lnx - 3ivr 



ff = [ l dxg { l\x)(?>- — — hix-Mii), (16) 
Jo V 1 — X J 

where fj has the same expression as the hard scattering kernel fx in B — > tttt ||. The 
hard kernel fj T for hard spectator interactions have the expressions (Vi: recoiled meson, V^- 
emitted meson): 

N c h Jo p Jo rj Jo £ 

/f±f = _|!M£ 2(lTl) f ip m f dn *JM f\ g ™ {() 

JM C mBti± Jo p Jo r] z Jo 

4n 2 2f B f Vimi /-i , $f(p) f 1 ,_,A VxMf-s V2lv),^Z + fj 



f dp^r dfjdd 9 ^\fj) g ^\o 

Jo p Jo 



N c m? B h± Jo ~' r p Jo '"" > Y'- L ^ 2 

,1 Vi(«)/-\ V 2 (a)/ f N £ + V _,_ 1 Vi(a)/-\ V 2 (v) / £\ + V \ 



where 
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h = (ml -m\- m 2 2 )(m B + m 1 )A^{w 2 2 ) - JllM^ A^(mj), 

TUB + Tn\ 

h ± = (m B + m^A^iml) T 2m f Pc V BVl (m 2 2 ), (18) 

and we have neglected the light quark masses and applied the approximation p ~ 0, and the 
B meson wave function 01: 



ifr 

4 iyr ' lu,H 'Jo 

with = (1, 0, 0, —1). Note that the presence of logarithmic and linear infrared divergences 
in ffj implies that the spectator interaction is dominated by soft gluon exchanges in the 
final states. Hence, factorization breaks down at the twist-3 order for transversely polarized 
vector meson states. We will introduce a cutoff of order AQCT>/?nb to regulate the linear and 
logarithmic divergences. The choice of the cutoff is not important here since the transversely 
polarized amplitudes are suppressed anyway. 

Two remarks are in order, (i) Since (V^^lO) = 0, B — > VV decays do not receive fac- 
torizable contributions from a 6 and a§ penguin terms except for spacelike penguin diagrams 
0. (ii) The first two terms — 12 ln(/i/mt,) — 18 in Eq. ( |l~5|) for helcity ±± states arise from 
the twist-3 DA g±\u) and will render the parameters af (except for a§ and a§ ) scale and 
scheme independent. 



4. To proceed we compute the branching ratios using LCSR and BSW models for heavy- 
light form factors (see Table I). The factorized amplitudes of B — > VV modes are given in 
IHII. Note that the original BSW model assumes a monopole behavior for all the form factors. 
This is not consistent with heavy quark symmetry for heavy-to-heavy transition. Therefore, 
we will employ the BSW model for the heavy-to-light form factors at zero momentum transfer 
but take a different ansatz for their q 2 dependence, namely a dipole dependence for A ,A 2 
and V. In the light-cone sum rule analysis, the form-factor q 2 dependence is given in M. 



To proceed we use the next-to-leading Wilson coefficients in the NDR scheme [12 



ci = 1.082, c 2 = -0.185, c 3 = 0.014, c 4 = -0.035, c 5 = 0.009, c 6 = -0.041, 

c 7 / a = -0.002, c 8 /a = 0.054, c 9 /a = -1.292, c 10 /a = 0.263, c g = -0.143, (20) 



with a being an electromagnetic fine-structure coupling constant. For the LCDAs, we use 



the asymptotic form for the vector meson |10 



gi\x) = 6 -[l + (2x-l) 2 



(21) 



and the B meson wave function 
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$f(p) = N B p\l-p) 2 exp 

with uj b = 0.25 GeV and N B being a normalization constant. For the decay constants, we 
use 

f p = 216 MeV, f K * = 221 MeV, f u = 195 MeV, / = 237 MeV, (23) 
and we will assume fy = jy for the tensor decay constant. 

TABLE I. Form factors at zero momentum transfer for B — > P and B — > V transitions evaluated 
in the light-cone sum rule (LCSR) analysis [3]. The values given in the square brackets are obtained 
in the BSW model [2]. We have assumed SU(3) symmetry for the B — > uj form factors in the LCSR 
approach. 



Decay 


V 


A x 


A 2 


^3 = A 


B -> p ± 


0.338 [0.329] 


0.261 [0.283] 


0.223 [0.283] 


0.372 [0.281] 


B -> u 


0.239 [0.232] 


0.185 [0.199] 


0.158 [0.199] 


0.263 [0.198] 


B -> K* 


0.458 [0.369] 


0.337 [0.328] 


0.283 [0.331] 


0.470 [0.321] 



1 r pm B Y 

2 V uj b J 



(22) 



To illustrate the non-universality of nonfactorizable effects for helicity amplitudes, we 
give a few numerical results for the parameters a\: 

o? = 1.04 + 0.01z, at = 1.02 + 0.0H, Oj = 1.11 + 0.04i, 

a° 2 = 0.09 - 0.08i, a+ = 0.17-0.08z, a" = -0.38 - 0.25z, 

a° = -0.033 - 0.004i, a+ = -0.026 - 0.004i, a7 = -0.040 - 0.007i, (24) 

in the LCSR model for form factors. Therefore, nonfactorizable corrections to helicity ampli- 
tudes are not universal. From Table II we see that the branching ratios predicted by LCSR 
is larger than that by the BSW model by a factor of 1.5 ~ 2. Evidently, the experimental 
results for B — > <j)K* favor the LCSR form factors for B — V transition. It should be stressed 
that thus far we have not taken into account power corrections such as annihilation diagrams 
and higher-twist wave functions for the longitudinally polarized vector meson. In particular, 
weak annihilations induced by the (S — P)(S + P) penguin operators are no longer subject 
to helicity suppression and hence can be sizable (see the last term in Eq. ([?]) and [0). How- 



ever, contrary to the PP and PV modes, the annihilation amplitude in the VV case does 
not gain a chiral enhancement of order m B /(m q mb). Therefore, it is truly power suppressed 
in the heavy quark limit. 



8 



It is also clear from Table II that the tree-dominated modes p + p~, p~p°, p~oo have 
larger branching ratios of order (2 ~ 3) x 10 -5 than the penguin-dominated ones. This 
is ascribed to the fact that the (S — P)(S + P) penguin operators do not contribute to 
factorizable W^-emission amplitudes. By contrast, the p°p° and ujuj modes have rather small 
branching ratios because the parameter 02 is small in QCD factorization. We have also 

computed |i?oo| 2 an d |-^±±| 2 for each channel and found that \H / Hoo\ 2 = (5 ~ 20)% and 

\H ++ /H 00 \ 2 = (10~ 5 ~ 10- 3 ). 

5. We have analyzed B — > VV decays within the framework of QCD factorization. 
Contrary to phenomenological generalized factorization, nonfactorizable corrections to each 
partial-wave or helicity amplitude are not the same; the effective parameters a* vary for dif- 
ferent helicity amplitudes. The leading-twist nonfactorizable corrections to the transversely 
polarized amplitudes vanish in the chiral limit and hence it is necessary to take into account 
twist-3 distribution amplitudes of the vector meson in order to have renormalization scale 
and scheme independent predictions. Branching ratios of B — > VV decays are calculated 
in two different models for form factors, and the predicted decay rates are different by a 
factor of 1.5 ~ 2. In the heavy quark limit, both vector mesons in the charmless B — > VV 
decay should have zero helicity and the corresponding amplitude is proportional to the form 
factor difference (A\ — A2). The recent observation of B — * <pK* indicates that the form 
factors Ai(q 2 ) and A^{q 2 ) cannot be very similar at low q 2 as implied by the BSW model. 
Owing to the absence of (S — P)(S + P) penguin operator contributions to W^-emission 
amplitudes, tree-dominated B — > VV decays tend to have larger branching ratios than the 
penguin-dominated ones. 
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TABLE II. Branching ratios (in units of 10 -6 ) averaged over CP-conjugate modes for charmless 
B — > VV decays. Two different form-factor models, the LCSR and the BSW models, are adopted 
and the unitarity angle 7 = 60° is employed. Experimental limits and results are taken from 
[14-16,4-6 and the limits indicated by * are quoted from [14] only for the helicity zero states. 



Decay LCSR BSW Expt. 

B U p~p+ 3ol) 2L2 < 2200 

B° -> p°p° 0.26 0.20 < 5.9* 

B°^ujuj 0.30 0.22 < 19 

B~ -> p~p° 21.8 12.8 < 120 

B- -» p-u 21.0 13.8 < 47 

B° -> K*-p+ 4.84 3.12 

B° -> 0.99 0.71 < 19* 

_> K*°K*° 0.32 0.16 < 10* 

B~ -» if*"/ 5.59 2.94 < 54* 

B~ -» 7T V 6.70 4.01 

B~ -> K*-R*° 0.34 0.15 < 50* 

B° -» p°(j) 0.003 0.002 < 13 

B°^uj(j) 0.003 0.002 < 21 

B- -> p-0 0.007 0.004 < 16 

B° ^ p°lo 0.12 0.07 < 11 

7J° _> k*°u; 3.66 2.15 < 19 

5- -> K*-^ 3.12 1.88 < 52 

B~ -» K*-(p 9.30 4.32 9.7±|;| ± 1.7 



B° -> K*V 8.71 4.62 8.6j|;| ± 1.1 



10.6 



+6.4+1.8 
-4.9-1.6 



< 36 



11.5 

;+8 



+4.5+1.8 
-3.7-1.7 



1511 ± 3 |5[ 
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